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The state of a non-relativistic gravitational dynamical system is known at any time t if the dynamical rule, i.e. 
Newton's equations of motion, can be solved; this requires specification of the gravitational potential. The 
evolution of a bunch of phase space coordinates w is deterministic, though generally non-linear. We discuss 
the novel Bayesian non-parametric algorithm CHASSIS that gives phase space pdf /(w) and potential <E>(x) of a 
relaxed gravitational system. CHASSIS is undemanding in terms of input requirements in that it is viable given 
incomplete, single-component velocity information of system members. Here x is the 3-D spatial coordinate and 
w = x + v where v is the 3-D velocity vector. CHASSIS works with a 2-integral / = f(E, L) where energy 
£■ = <!> + v 2 /2, v 2 = 2? =1 v? and the angular momentum is L = |r x v|, where r is the spherical spatial vector. Also, 
we assume spherical symmetry. CHASSIS obtains the /(•) from which the kinematic data is most likely to have 
been drawn, in the best choice for <!>(■), using an MCMC optimiser (Metropolis-Hastings). The likelihood function 
X is defined in terms of the projections of /(■) into the space of observables and the maximum in X is sought 
by the optimiser. The recovered solutions can be susceptible to large uncertainties given the dimensionality of 
the domain of the unconstrained /(■) and the typically small, observed velocity samples in distant astrophysical 
systems. This scenario is tackled by assuming / = f(E), i.e. we assume the phase space to be isotropic. 
However, this simplifying assumption of isotropy is addressed by undertaking a Bayesian test of significance that 
is developed to be used in the non-parametric context. 

A test based on the /j-value estimates of the goodness of isotropy in the data was previously undertaken 
(Chakrabarty & Saha 2001). However, p-values are sensitive to sample sizes and obfuscate interpretation of 
analyses of differently sized kinematic samples. Thus, a Bayesian formalism is a better alternative, eg. Fully 
Bayesian Significance Test or FBST (Pereira & Stern 1999, Pereira, Stern & Wechsler 2008). The null hypothesis 
that we aim to test, is that the data are drawn from an isotropic /(■), i.e. Hq : / = v F[£ , (v 2 /2 + <I'(r))] where the data 
are drawn from / and *P is some function: > for E < and ^ = otherwise. Within FBST, the evidence value 
(ev) in favour of Hq is obtained by first numerically spotting the most likely configuration (6*) that is compatible 
with //(> and then finding the volume of the tangential set T by numerical integration. Here T is the set of all 
configurations with posterior probability in excess of that of 6*. We have developed the implementation of this 
scheme in the non-parametric context; in CHASSIS, the configurations are f(E) and <£>(r). To have configurations 
obeying Hq, we perform sampling from the f(E) - <D(r) pair identified upon convergence of a run of CHASSIS. 
From this sampling, the resulting f(E) - O(r) configuration corresponding to the highest _£ is compared to all the 
other f(E) - O(r) pairs, in order to obtain a measure of the volume of T. 

We discuss two distinct applications of the isotropic version of CHASSIS. In one, 2 distinct kinematic data 
sets of 2 distinct types of members of an example galaxy are analysed by CHASSIS under the assumption of 
isotropy. The /(•) and <£(•) recovered from runs done with the two data sets are identified as distinct. Given that 
the same galaxy cannot be described by two different gravitational potentials, the risk involved in the very method 
of extracting the galactic potential from kinematic data of individual galactic members is demonstrated here for 
the first time. The goodness of the assumption of isotropy, given the 2 data sets, is quantified by our Bayesian 
test of hypothesis. 

In the second application, it is shown that once the amount of gravitational matter inside a fiduciary radius 
is pinned down from independent measurements, the recovered <£(•) is unique, irrespective of the toy /(•), from 
which kinematic data samples are drawn as input for CHASSIS, as long as velocity dispersion values are mea- 
sured at 1 or more different locations in the system. This consistent nature of (£>(•) is arrived at, notwithstanding 
varied forms for the assumed toy /(•), including isotropic as well as E & L dependent forms. 

Keywords: Bayes theorem, Bayesian significance test, Astrophysical applications. 
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1 Introduction 

The complete characterisation of a gravitationally bound, non-relativistic dynamical system can be undertaken 
with the help of /(w, f) - the pdf of phase space W - and the gravitational potential <£(x); here t is time and 
w = x + v, where x represents the spatial coordinates and the velocity vector is v = x. A sample of phase space 
coordinates can be drawn from /(•) and allowed to evolve in <£(•), in accordance to Newton's laws. In this way, 
the evolution of the system is deterministic at any time t, though non-linear in general. Hence, we aim to estimate 
/(•) and <t>(-); we focus on the characterisation of astrophysical systems in this paper. A related aim is to derive 
the distribution of the total gravitating matter from the estimated <&(■), while keeping in mind that such total mass 
is accounted for only partly by luminous matter while the greater fraction is dark matter in these systems. 
Conventionally, 

• /(•) and <£(■) are almost always parametrically described. However, given that astrophysical systems such 
as galaxies are more likely to manifest complexity in their dynamics than otherwise, any smooth parametric 
description of such systems is erroneous. 

• mass determination is typically pursued via observed photometric or luminous information though no func- 
tional dependence of the total (luminous +dark) matter content on such measurements exist. 

• inhomogeneities in the measurement errors notwithstanding, goodness-of-fit parameters are often invoked 
to seek the solution. For galaxies, measurements are typically noisy and such goodness-of-fit parameters 
can be artificially inflated (Bissantz & Munk, 2001). 

All these issues suggest a better - preferably, a nonparametric - route to /(•) and <&(•) determination. It is such a 
novel, data driven characterisation of real (as distinguished from simulated) astrophysical systems that CHAS- 
SIS offers, using the few kinematic measurements that are typically available. Importantly, we test the chief 
assumption of our algorithm using a test of significance that is developed in this regard. 

As motivated above, we discard all photometric information that may be available for a system at hand and 
use the kinematic information that is sometimes available, namely velocities along the line-of-sight (LOS) of 
individual system members, (such as stars). We define the galaxy such that the x\ - X2 coordinate system spans 
the plane-of-the-sky (POS) and the LOS is along the z-axis. Thus, our data comprise 1 -component V3 values of 
individual galactic members and their POS coordinates. We also account for the errors in V3 measurement. The 
observed data samples often bear < 100 data points. These data are input into the Bayesian algorithm CHASSIS 
(Chakrabarty & Saha 2001, Chakrabarty & Portegies Zwart 2005). 

2 CHASSIS 

CHASSIS helps constrain <&(•) and /(■) of galaxies, given the data u (say). Actually, within CHASSIS, we 
constrain the gravitational matter density p(-) rather than <£(■), where Poisson equation connects p(-) and <!>(■) as 
in: 

O(x) = -4ttGV 2 p(x) (1) 

This helps to avoid problems about negative p(-). The calculation of O(-) from p(-) is undertaken at every step. 
Dynamical theory tells us (Dubrovin, Fomenko & Petrovich, 1990): 

/ = f[Ki(w)] , where K t = Vz = 1 , 2, 3 . . . . (2) 

i.e. Kj is an integral of motion. Now, we realise that the size of the data u is small to moderate and typically 
bears no more information other than a single component of velocity. In such a case, we feel that the data are 
not sufficient to constrain the extended forms of /(•) and <!>(•)• In other words, we resort to making assumptions 
about /(•) and <&(•). 

In fact, we assume Kj ^constant, only for i = 1,2. Thus, / = f[K[{w), ^(w)] where K\ =energy E of the 
galactic particle and K2 = the angular momentum or L. Here E = J^j vj/2 + €>(x) and L = |r x v] where r is the 
spherical radius and v is the 3-D velocity vector: J^j x 2 = r 2 , v 2 = Ytj tf- We also assume radial symmetry in the 
potential, i.e. 4>(x) = O(r). Thus, the particle energy is E = v 2 /2 + O(r). In this geometry, the Poisson equation 
(EquationQ]) is solved numerically by assuming the mass to be stratified on spherical shells. In other words, the 
relevant radial range is discretised and over each radial bin, p(r) is held a constant. 

In this background, we seek 

Pr(/,p| M )(xPr( M |/,p)Pr(/,p) (3) 
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where the only constraints posed by the priors are physically realistic requirements of positivity and mono tonicity. 
This owes to the fact that in general, we do not possess any other prior information that could help constrain the 
forms of the sought functions. Thus, 



Pr(/,p) = Pr(/)Pr(p), where 

< 

an L 

otherwise 



Pr(/) = 1 if / > V % 

oE 



Pr(p) = 1 if p > V — < 

dr 

= otherwise (4) 

We adopt the above priors and estimate the /(•) from which the data u is most likely to have been drawn, in 
the estimated <!>(•)■ This is done by iterating towards the most likely set of {/(■), p(-)} starting with an arbitrarily 
chosen seed. At every step, the current choice of /(•) is projected into the space of observables (spanned by 
xi,x 2 , V3); a ready definition for the likelihood function X is in terms of the projection r](-) of f(E, L). 

r]i(x\,x' 2 ,v\) = J f[E(v, <£(>)), L] dx 3 dv l dv 2 and 

Ndata 

L = £ log 77, (5) 

i=i 

where x\, x' 2 , v' 3 is the i th data point in the Ndata sized data sample. 

The numerical implementation of a trial /(•) function, at a trial (p(-) or) <!>(•) is the crucially important question 
from the point of view of algorithm design. We do this by discretising the E-L space and holding /(•) a constant 
(=fceii) over a given E-L cell. The contribution to 77O), from this E-L cell - defined by say, E e [E\,E 2 ], L e 
[L\,Li\ - is given as: 



n? 1 = Uu J 'dx 3 dndv 2 and 



cell 

Let the integral on the right hand side of the former of these two equations be A. Then, we seek the mapping 
A : — > E - Lspace. 

In order to establish this, first we determine the 2-D area of intersection between the locii of E — E\, E — E 2 , 
L - L\, L — L 2 , in the v\ - v 2 space. This gives the connection between the E-L cell at hand and v\,v 2 . 
Mapping X3 to this cell requires knowledge of the minimum and maximum values of X3 that are allowed in 

the cell, given the data point x\, x' 2 , v\. This maximum value is - (x'j) 2 - (x l 2 ) 2 where ro is the solution to: 

E 2 = v\/2 + L 2 /2rl + fl>(r ). The minimum x 3 is 0. 

This explains the background to the structures of p(-) and /(•). 



2.1 p-histogram and /-histogram 

The representation of f(E, L) over the discretised E-L space, is akin to a 2-D histogram. Similarly, the p(r) 
structure is represented as a 1-D histogram. These histograms are updated at the beginning of every step, while 
maintaining positivity and monotonicity. 

The jump distribution we use is discussed below. If in step k, for r e [r q -i,r q ] (Vg = 1, . . . ,N r , ro=0), 
p(r q ) = a k q , then in the k + I th step: 

exp(^) (7) 

where % is a random number with % € [-0.5,0.5] and s\ is an experimentally optimised scale that determines 
the scale over which the shape of the p-histogram is changed. This updating is done Vq. Once the shape of the 
p-histogram is updated in this way, the whole structure is scaled by the factor e\pCR/s 2 ) where 52 is another scale. 
The /-histogram is similarly updated in shape and subsequently normalised. 
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2.2 Optimisation 

Once the histograms are updated, we project the current /(•) over an E - L cell, into the space of observables, 
for the i' h data point and then sum over all such cells to get m (Equation 7). This is done V i, to obtain £ 
(Equation 6). The global maxima in £ is sought by the Metropolis-Hastings algorithm (Metropolis et. al 1953, 
Hastings 1970, Chib & Greenberg 1995). Anticipating the likelihood distribution to be multimodal, we work with 
highly dispersed seeds to initiate several chains (Gelman & Rubin 1992) as well as employ simulated annealing 
on a single chain. The latter approach, though perhaps less obvious, is one that we find very effective in test runs. 

While the optimisation routine is hard-wired within CHASSIS, the user is allowed the flexibility to adjust 
details of the used cooling schedule and other optimisation parameters such as the scales relevant to the jump 
distribution (s\, s 2 )- In this note, it is worth mentioning that the current implementation of optimisation is 
modular, and it is simple to replace it by a more effcient routine. 

2.3 Required User Input 

The methodology discussed above is incorporated into CHASSIS and all that user is required to input is the 
velocity data, the source of which is independent of CHASSIS. Thus, measured kinematic data, irrespective of 
its source, is acceptable, as long as the columns pertain to the observables r p , V3 and the measured errors in V3. 

Here r p = ^jx 2 x + x 2 v Besides, the user is allowed to input details such as the number of bins, bin widths, fraction 
of data she wants to perform the run with, the seeds for the sought solutions and the optimisation related details 
(see Section 2.2). The user inputs are advanced via an input file that CHASSIS calls at the beginning of a run. 

2.4 Assumption of Isotropy 

Given the limited data sample, we find that limiting the domain of / to 2-D is not constraining enough in reducing 
the magnitude of uncertainties in the estimated solutions to useful levels. Thus, we resort to imposing the further 
constraint that / = f(E), i.e. we assume isotropy to exist in phase space (since E is symmetrical in vj and xj, 
where 7=1,2,3). Then we resort to (1) justifying or rejecting our assumption in the data by performing a test of 
hypothesis exercise (2)exploring independent measurements that may be available in the literature to obtain a p 
that is unaffected by the amount of anisotropy in the data. 

3 Testing for Isotropy - Nonparametric FBST 

We test for isotropy in the data, i.e. the null hypothesis Ho : f = x P[£'(v, r)], where / is the phase space pdf from 
which the observed data are drawn and ¥ is some function that manifests phase space isotropy: *¥(E) = 0V£ > 
and > otherwise. 

This Hq is tested in the data u along the lines of the Fully Bayesian Significance Test or FBST (Pereira & 
Stern 1999; Pereira, Stern & Wechsler 2008), except that here, we advance a nonparametric implementation of 
the same. Our null is sharp, as is the requirement for FBST (Madurga, Esteves & Wechsler 2001). We refer to our 
object functions {p( •),/(•)} = 9 ( sav ); let 9 e 0-space. We assume that 9 is continuous in the ©-space. According 
to FBST, the evidence in favour of Ho is 1 - ev, where: 

ev = 1 - Pv(9 e T\u), where 

T = {9 : Pr(9\u) > Pr(9*\H )}. (8) 

Here 9* is the value of 9 which, under the null, maximises the posterior Pr(9\u). 

At the end of every iterative step during a run of CHASSIS, a 9 configuration is identified. The /(•) recovered 
upon convergence of the run is indeed a function of E and E only but the true / is not necessarily so. 

Upon convergence of a run of CHASSIS performed with data u, we sample the recovered 9, N times, such 
that the i th sampling of 9 gives the i th set of observables or uf, i = 1,2, ...,N. These N data samples are then 
input into N different new runs of the algorithm. During the run performed with the data w,-, the /* iterative step 
yields the configuration 9\ (say), where 9\ = {/(•); >P(");}- Then /(•)■ is isotropic V;', j, since the phase space pdf 
from which m, is drawn is the recovered f(E), which by construction, is indeed isotropic. 

We scan over all i, j to identify i*,f for which the posterior is maximised. Thus, 9j, = 9* is identified. Here 
6* = [f ""(E), p*(r)}, i.e. functions recovered at the end of the j*-th step, in a run performed with the data 

In this nonparametric implementation, Pr(6> e T\u) = X/Y, where X is the number of times that a step yields a 
likelihood in excess of Pr(0*|// o ) m a H the undertaken runs. Y is the total number of iterative steps in all the runs 
undertaken. 
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Figure 1: Left: The (normalised) phase space distribution functions recovered from two runs of CHASSIS, 
performed with the GC data uq (in black) and PNe data up (in red). The normalisation is performed to ensure 
that f(E)= 1 for E=-\. These pdfs are recovered under the assumption of isotropy and our implementation of 
nonparametric FBST shows that the profile in black is expected to be closer to the true phase space pdf than is the 
profile in red. The errors are ±l-<x uncertainties identified on the solution by the optimiser. Middle: Gravitational 
matter density distributions from three runs with u c performed with three different seeds (in magenta, cyan and 
black) and from 3 runs with u P (in yellow, red and blue). The units of density and radius are astrophysical. Right: 
p(r) from one of the runs done with «c is shown in black while that with up is in red. The p*(r) corresponding to 
the implementation of data uq is shown in green while that corresponding to up is in blue. We notice that pc(r) 
is consistent with p*(r) obtained from uq while p*(r) obtained from up is significantly lower than pp(r). 

4 Application 1 - multistability in galaxies 

It is a common practise in astrophysics to employ the measured 1-D velocity data of suitable galactic members, 
with the aim of recovering the total gravitational mass distributions. To examine the viability of such a practise, 
we employ the available measured kinematic data of two distinct classes of galactic members in an example 
galaxy - as an aside, these are planetary nebulae (PNe) and globular clusters (GCs). The data of 164 PNe (u P ) are 
due to Douglas et al (2007) and that of 30 GCs (u G ) are due to Bergond et. al (2006). These sample sizes are too 
discrepant to allow for easy interpretation of any p-value based testing of Ho defined in Section 3. Instead, we 
resort to the non-parametric FBST discussed above. 

The two data sets are input into the isotropy-assuming and sphericity-assuming CHASSIS. The f{E) and p(r) 
recovered from the two distinct data are found to be inconsistent with each other, within error bars (Figure 1). 
The difference in the recovered /(•) could arise from different divergences between an isotropic phase space pdf 
and the / from which up are drawn, as compared to that from which uq are drawn. However, the distribution of 
gravitational matter in the galaxy should be uniquely determined. That such is not our conclusion, prompts us to 
examine if the inherent assumption of isotropy is to be blamed. Thus, we test for our null Hq (defined above in 
Section 3) in the data up and uq separately. 

The results of our implementation of FBST are shown in Figure 1. We find that for three different runs done 
with distinct seeds, given u P , the average evidence in favour of the null is about 0.60 For three runs done with 
different seeds, given uq, the average 1 — evis 0.95. Thus we conclude that the degree of isotropy of the pdf that 
uq is drawn from is higher than that of the pdf that up is drawn from. 

We wonder if the difference in the recovered mass distributions be due to the concluded difference in isotropy 
in the two data sets? To understand this, we invoke the peculiarity of CHASSIS that the algorithm overestimates 
mass density at all r where phase space anisotropy prevails (discussed in Section 5). Thus, we expect that p(r) 
recovered using up (pp(r)) is more of an over-estimate compared to the galactic mass density than is the p(r) 
recovered using uq {pair)). 

However, as shown in Figure 1, at all r > 6 kpc, pc(r) > pp(r). Therefore, to reconcile the difference between 
pp(r) and pc(r), the isotropy issue cannot help unless we propose that uq is drawn from a more anisotropic pdf 
than up. This is of course not true but its inverse is. Hence we conclude that differences in pc(r) and pp(r) are 
intrinsic to the system and not due to our assumption of isotropy. 

Thus we have demonstrated the potential risk in employing kinematic data of individual galactic members of 
a particular population type, to compute the gravitational mass distribution of galaxies. We have also shown that 
the galactic phase space is described by at least two distinct basins of attractions, i.e. the galaxy is multistable, as 
we would expect complex systems like galaxies to be. 
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Figure 2: Gravitational mass density distributions over radius, obtained from runs done with data drawn from 
three different phase space distributions (fcauss, Iwd fMkhw) that are distinguished from each other in terms of 
the inherent degree of anisotropy in their forms. p(r) estimated from run performed with data S i SO tropy is m black, 
with Swd is in red and with S Michie is in blue. The p(r) implied by the fiduciary test potential O teJ/ (r) is in cyan. 
The left panel displays results obtained when the constraint of total mass is included (Mq as 4.06 x 10 11 M within 
about 8.7 kpc). On the right, the general inconsistency between the estimated profiles, when the mass constraint 
is excluded, is brought out. 



5 Application 2 - using the total mass constraint 



Motivated by the need to simplify our analysis by reducing the number of degrees of freedom to the bare mini- 
mum, we persist with the assumption of / = f{E), i.e. phase space is isotropic. We envisage that when the data 
have been drawn from an anisotropic pdf, the algorithm will imply erroneous answers. Using physical arguments 
we can predict the nature of this error - CHASSIS overestimates p(r) at r where phase space anisotropy prevails. 

We implement an independent measure of the total gravitational matter (Mo) inside a given radius (Re) in 
an example galaxy with the aim of recovering the correct solution for p(r) under our assumption of isotropy, 
irrespective of the degree of anisotropy of the pdf from which the input data is drawn. If the system is at a large 
distance from us, we cannot get velocity data of individual galactic members. Then, we can only get projected 
velocity dispersion values (cr p ) at 1 or a few radial locations in the galaxy. In this background of sparse and 
incomplete velocity data, we prepare velocity data samples for inputting to CHASSIS, in the following way. 

We select observables (x\,xi, V3) from 3 different toy /(•), two of which are selected to depend on E and L 
while the third is isotropic. For the example galaxy, cr p is known at r\, r2 and say. Then we consider the system 
to be divided into 3 anulii with r e [0, r{\, (r\,r-{\, (ji, r^]. The error 6 in the cr p measured at r,- = 1, 2, 3) can 
be related to the size of the data sample Nj we intend to draw from the given annulus, assuming normal error 
distribution. The phase space pdf that we choose our samples from are fcauss, fwD fuichie- 



fcauss(E) — 



fwD(E, L) 



1 



6) 



fMichie(E, L) — 



-exp 
lira 2 
E >0, 
1 

— exp 

E >0. 

1 / L 2 

exp ; 

V27TO- 2 \ r a <r 



E < 0, 



exp 



E < 0, 



E < 0, 



= E > 0. 



(9) 



The samples chosen from fcauss, fwD fuichie are Scauss, S wdS MicHe (say). Also, to test for the effect of data from 
different forms of /(•), we define E = v 2 /2 + O fes ,( r ) where we choose Q>test(f) ~ 1/ -y/r 2 + r 2 . 

The total mass constraint is expected to narrow down the range of solutions possible; in this sense it acts as 
a prior on the solution for p(-). We test if the chosen data samples, when input into CHASSIS, recover a cD(r) 
that is concurrent with (£> tes ,(r) when we include/exclude the constraint that total mass within r — Re is Mq. 
Here M is obtained from literature as about 4.06xlO n M , with errors of SM = +0.2 x 10 11 M Q and R E -8.7 
kpc (Koopmans & Treu 2003). We incorporate the constraint by adding a penalty function to the definition 
of the likelihood; the role of this penalty function is to penalise any solution that implies a total gravitational 
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mass within R E (M C (R E )) different from Mo. This penalty function is a\M c (R E ) - Mq\/26Mq. Here a is a flag, 
designed to include or exclude the constraint from the definition of the likelihood, depending on whether it is 1 
or respectively. 

The results of conducting our test runs are shown in Figure 2. We find that when the constraint is included, 
the p(r) profiles are consistent with each other within errors, irrespective of the anisotropy in the data used to 
obtain this profile. In absence of guidance from the constraint, anisotropy affects results. 

6 Summary 

Here we have discussed the novel nonparametric Bayesian algorithm CHASSIS that estimates the most likely 
phase space pdf (/(•)) from which an observed sample of 1-D component of velocities of individual galactic 
members is drawn, at the most likely gravitational potential <!>(•) of the galaxy. The main purpose of this paper 
is to bring out the fact that CHASSIS is a robust and viable algorithm that can be implemented to extract the all- 
important gravitational matter density distribution in distant galaxies, even within the domain of very sparse and 
incomplete data. Given the dearth of measurements in these systems, it is prudent to work with a small number 
of degrees of freedom. With this in mind, one version of CHASSIS has been designed to work under the purview 
of the assumption of isotropy in phase space, by which we imply an /(•) that is symmetric in the 3 velocity and 
3 spatial coordinates. In another version, the assumption of isotropy is not made and CHASSIS is made to work 
with a greater number of dof (Chakrabarty & Saha, under preparation). 

We offer independent means of tackling the obvious fallout of the assumption of isotropy, when invoked; 
phase space pdfs from which realistic data are drawn, will not be isotropic in general, leading to spurious solu- 
tions for /(•) and <!>(•)• Firstly, a robust test of hypothesis is developed that tests the assumption of isotropy, given 
the data. This test is modelled after Pereira & Stern's (1999) FBST and is designed to work in the nonparametric 
context. Data available in astrophysical literature are implemented to conclude that kinematic data of distinct 
galactic member populations will in general provide distinct gravitational mass distributions. The multistability 
of the example galaxy is demonstrated and the folly of this mode of mass determination is indicated. Secondly, 
we demonstrate that the usage of information about total gravitational mass, available in the literature, can con- 
strain the sought gravitational matter density distribution, irrespective of anisotropy in the data. Such a constraint 
basically supplements for the uninformative priors that we use within CHASSIS. 

Currently, we are working on the establishment of a critical value of the evidence value in favour of a given 
null. This will enable the quantified judgement of when to reject or accept the null, given the data. At the moment, 
our implementation of FBST only allows for a comparative judgement. 
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